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Abstract 

Using the CLEO II.V detector observing e + e~ collisions at around 10.6 GeV we search for neutral 
D mixing in semileptonic D° decays tagged in charged D* decays. Combining the results from the 
Keu and K*eu channels we find that the rate for D mixing is less than 0.0078 at 90% C.L. 
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The study of mixing in the K° and sectors has provided a wealth of information to 
guide the form and content of the Standard Model. In the framework of the Standard Model, 
mixing in the charm meson sector is predicted to be small |l( , making this an excellent place 
to search for non-Standard Model effects. 

A D° can evolve into a D° through well known, on-shell intermediate states, or through 
off-shell intermediate states such as those that might be present due to new physics. We 
denote the amplitude through the former (latter) states by — iy (x), in units of T D o/2 @. 
The Standard Model contributions to x are suppressed to \x\ ~ tan 2 #c* ~ 5% and the 
Glashow-Illiopolous-Maiani ||| cancellation could further suppress \x\ down to 10~ 6 — 1CT 2 . 
Many non-Standard Model processes could lead to \x\ > 1%. Signatures of new physics 
include \x\ 3> \y\ and CP violating interference between x and y or between x and a direct 
decay amplitude. 

Observation of D mixing in hadronic decay channels is complicated by doubly Cabibbo 
suppressed decays, where both the decay of the original charm quark and subsequent charged 
W decay proceed in Cabibbo suppressed modes. Such a decay mimics mixing of a D° to 
a D° followed by the dominant Cabibbo allowed channel decay. In semileptonic decays no 
such double suppression is allowed as the W decays to a charged lepton and neutrino, and 
the charge of the lepton tags whether a c, producing a positive lepton, or c, producing a 
negative lepton, has decayed. When the production flavor of the D is tagged in charged D* 
decay or some other way, a wrong sign lepton produced in a subsequent semileptonic decay 
unambiguously signals D mixing. Other mechanisms that produce leptons of the opposite 
sign in D° or D° decay are highly suppressed. The integrated mixing rate normalized to 
the total decay rate is equal to |(x 2 + y 2 ), and is called Rm- This is at once good and bad 
news. The observation of any wrong sign semileptonic D decay is an unambiguous signal 
of D mixing. With both x and y small (< (9(0.01)) the rate of wrong sign semileptonic D 
decays will be very small (< (9(0.0001)), and an observation of D mixing in semileptonic 
decays would not give insight on the relative sizes of x and y. In contrast hadronic decays 
have wrong sign contributions from both doubly Cabibbo suppressed decays and mixing. 
The two channels interfere resulting in a term that depends on a linear combination of x 
and y, and the proper time dependence of wrong sign final states can be used to measure x 
and y separately. 

The proper decay time dependence of semileptonic mixed final states in units of the mean 
D° lifetime, t D o = (410.3±1.5)fs 0], is r(t) = (x 2 +y 2 )t 2 e~ t . Thus mixed semileptonic decays 
should populate larger proper times with an average decay time of three. The dominant 
direct decays are distributed as e _< with an average decay time of one. This difference is 
taken advantage of to increase the sensitivity to D mixing. 

To date no one has observed evidence for D mixing in any channel. The best limit on 
Rm comes from the FOCUS collaboration [^| where they find Rm < 0.0010 at 90% C.L. 
Belle Q reports a similar limit, Rm < 0.0014 at 90% C.L. and BABAR 0] a higher limit of 
Rm < 0.0042 at 90% C.L. Other limits on D mixing are summarized in Our experimental 
situation is very similar to Belle and BABAR, but our data sample is more than an order 
of magnitude smaller. Thus we expect not to be sensitive to any D mixing signal and set 
limits higher than above. 

Our data were collected using the CLEO II. V upgrade [j| of the CLEO II detector [l(j 
between February 1996 and February 1999 at the Cornell Electron Storage Ring (CESR). 
The data correspond to 9.0 fb _1 of e + e~ collisions near a/s « 10.6 GeV. The detector 
consisted of cylindrical tracking chambers and an electro-magnetic calorimeter immersed 
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in a 1.5 Tesla axial magnetic field, surrounded by muon chambers. The reconstruction of 
displaced vertices from charm decays is made possible by the addition of a silicon vertex 
detector (SVX) in CLEO II. V. The charged particle trajectories are fit using a Kalman filter 
technique [ll[ that takes into account energy loss as the particles pass through the material 
of the beam pipe and detector. Specific ionization for charged particle identification is 
measured in the main drift chamber with a resolution of about 7%. Electrons above 500 
MeV/c momentum are also identified by matching of track momentum with calorimeter 
energy deposition and requiring that the calorimeter shower has the characteristics expected 
of an electro-magnetic rather than a hadronic shower. Hadrons are misidentified as electrons 
at roughly the 0.1% level from studies of known hadron samples in the data. Muons are 
not used as the CLEO muon identification system cannot cleanly identify muons below 
about 1.5 GeV/c momentum and the small sample of clean muons are not useful to study 
semileptonic D decays at our beam energy. 

To study backgrounds and relative selection efficiencies we use a GEANT El based de- 
tector simulation of our data. We use simulations of e + e _ — > qq with the quarks fragmenting 
and particles decaying generically guided by previous measurements. The data of this generic 
simulation corresponds to roughly ten times the luminosity collected by the detector. We 
also use simulations of e + e~ — > t + t~ and e + e~ — > BB to study small contributions to the 
background. For these studies the simulated samples are reconstructed and selected using 
the same methods as the data sample as described below. 

To select hadronic events and ensure that the event production point is well known there 
must be at least five well reconstructed tracks consistent with coming from the interaction 
region. The tracks must carry more than 15% of the collision energy. This selection is nearly 
100% efficient, removes events that would not pass subsequent reconstruction requirements, 
and ensures a well measured D° flight distance. All tracks used in the reconstruction of the 
decay chain except the pion from the charged D* decay are required to hit at least two of 
the three SVX planes in projections both transverse and parallel to the beam direction. 

The D°s are tagged at production in the decay of the charged D* to a charged pion and 
a D°. A 7r + indicates that a D° has been produced and a 7r~, a D°. Subsequent direct 
semileptonic decays of the D° or D° produce a charged lepton of the same sign as the pion 
from the D* decay, called right sign (RS) combinations, while semileptonic D decay after 
mixing produces a charged lepton of the opposite sign, called wrong sign (WS) combinations. 
Since little energy is released in the D* decay the pion is limited to a momentum of 400 
MeV/c at our beam energy, and is called the soft pion. We choose combinations such that 
the momentum of the D° is larger than 2.0 GeV/c and thus are dominated by D°s produced 
in e + e~ — > cc with only a small contribution produced by B decays. The analyses find RS 
combinations maximizing the signal to noise for the decay chains D* — > n so { t D° — > K^ev. 
The same selections are then used to find WS combinations, and the ratio of WS to RS 
combinations after accounting for background is a measure of Rm- The analysis of the 
D° — > Kev channel is described in full detail in 0] while the D° — > K*ev channel is 
similarly described in . 



The D° — > Kev analysis uses a neural net to distinguish signal events from background. 
All combinations of electrons, kaons identified loosely via specific ionization, and 7r so f t are 
considered. Eighteen variables are inputs to the net. They are selected such that our 
simulation describes their distribution well, and in general they describe the kinematics of 
the 7T so f t , kaon, and electron candidates separating random combinations from those produced 
in the desired decay chain. For example one input is the cosine of the angle between the 
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FIG. 1: The output of the neural net used in the D° — ► Kev analysis for the RS (left) and WS 
(right) samples comparing data (squares) with simulation (full histogram). Also shown with a 
dashed histogram is the predicted contribution of the signal. The very small WS signal is at the 
central value found by this analysis. 



electron and the kaon which peaks near one for the signal and at both one and negative one 
for background. The net is trained on our simulation and produces an output near one for 
signal-like combinations and negative one for background-like combinations. Figure Q shows 
the output of the neural net on RS and WS candidates. According to the simulation the RS 
sample is roughly 42% signal while we expect the WS to be dominated by background. We 
compare the output of the neural net thoroughly with the prediction of the simulation by 
varying both its input and structure. The simulation is found to predict very well changes 
observed in the data. Combinations with a neural net output of greater than 0.95 are 
selected for further analysis. 

To measure the decay time of the D° we refit the tracks requiring that the kaon and elec- 
tron come from a common D° decay vertex, use the thrust axis of the event as the direction 
of the D°, and require that the D° and 7r^ ow or D° and vr^ ow , come from a common vertex 
constrained to be in luminous region. This procedure improves the resolution on the decay 
time by 30% to about half a D° lifetime. Figure 121 shows the distribution of the decay times 
for the RS and WS samples. Overlaid are fits to a signal plus backgrounds. Background 
shapes are taken from the simulation and are dominated by non-charm events with a sig- 
nificant fraction of non-signal charm decays. The RS distribution is dominated by signal, 
while the WS is consistent with no signal and thus is dominated by background. Signal 
shapes are also taken from the simulation. The agreement between data and simulation is 
good and the normalizations for the RS signal and the background in both the RS and WS 
distributions agree with the predictions of the simulation. 

These fits find 2840±300 RS signal and 31±21 WS signal combinations. In the D° — > Kev 
channel we measure Rm = 0.0110 ± 0.0076 and see no evidence for mixing. Systematic 
uncertainties are discussed below. 

The D° — > K*ev analysis uses a more traditional approach. The charged K* is recon- 
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Proper Time (fs) Proper Time (fs) 

FIG. 2: The decay time for the D — > Keu analysis for the RS (left) and WS (right) samples. 
The squares show the data, the full histogram shows the fits, and the dashed histogram the signal 
contribution to the fits. 

structed in the K s ir channel followed by the K s decay to two charged pions. These two 
pions are fit to a common vertex and the mass of the tttc must be within 16 MeV/c 2 of 
the expected K s mass. Similarly for the K*, the K s -tc combination must be within 60 
MeV/c 2 of the known mass. When the K* is combined with an electron, and the pair is 
consistent with coming from a D° decay, the sample is fairly clean. The direction of the D° 
is determined with a weighted average of the thrust, the 7r so ft, and i^*-electron combination 
directions. This gives a better estimate of the energy release, Q, in the D* decay and sig- 
nal is distinguished from background in a fit to the two dimensional Q versus decay time 
distribution. The decay time is improved by a refit similar to the one described above. 

Figures E2 and 0] show the projections of the two dimensional fits on the Q and decay 
time axes for RS and WS candidates respectively. Shapes for signal and background are 
taken from the simulation. Agreement between the data and simulation is good both in the 
signal dominated RS sample and the background dominated WS sample. The prediction of 
the simulation is checked using D* — > n s \ ow D — > K s mnT° decays found in the data. The 7r° 
is ignored and the same methods are applied as in the K*ev analysis to reconstruct a two 
dimensional Q and decay time distribution. The simulation and the data agree very well in 
this check sample. The simulation also accurately predicts the size of the observed signal 
and background in both the RS and WS samples. 

These fits find 638±51 RS signal and — 30±8 WS signal combinations. When constrained 
to find at least zero WS signal the fit returns 0±2 WS signal combinations. The two results 
yield similar upper limits on Rm, and the latter is used to combine with the Kev channel 
as it yields a slightly more conservative upper limit. In the D° — * K*ev channel we measure 
Rm = 0.0000±0.0031 and see no evidence for mixing. Systematic uncertainties are discussed 
below. 

The two analyses are statistically independent and consistent; therefore they can be 
combined. We combine their central values weighting by their statistical uncertainties to 
find R M = 0.0016 ± 0.0029. 
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FIG. 3: The Q (top) and decay time (bottom) distributions for the D° — > K*ev analysis for the 
RS sample. The squares show the data, the histograms show projections of the two dimensional 
fit, and the dashed histogram shows the background contribution to the fit. 



The two analyses share some common systematic uncertainties. They both use the same 
simulation to model backgrounds and signal shapes. We take the smaller of the two sys- 
tematic uncertainties as the uncertainty on the combined value. The statistical uncertainty 
on simulated signal and background shapes causes an uncertainty of ±0.0023 on R M in the 
K*ev analysis and ±0.0028 in the Kev analysis. The shape of the background in the decay 
time is parametrized from the simulation in a similar way in both analyses. Variations in 
this parameterization affect the number of RS and WS signal combinations. This variation 
causes an uncertainty of ±0.0014 on Rm in the Kev analysis and ±0.0018 in the K*ev anal- 
ysis. The lower of the two uncertainties, the first in each case, is taken as the uncertainty 
on the combined result. 
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FIG. 4: The Q (top) and decay time (bottom) distributions for the D° — ► K*ev analysis for the WS 
sample. The squares shows the data and the histograms shows projections of the two dimensional 
fit. 



For systematic uncertainties not shared by the two analyses we add them by weighting 
their contribution to the central value. Specifically uncertainties in the Kev analysis con- 
tribute 15% of their size to the combined analysis, while uncertainties in the K*ev analysis 
contribute 85% of their size. In the K*eu analysis variations of the parametrization of the Q 
shape for signal and background add ±0.0008. In the Kev analysis the largest contribution 
to the uncertainty comes from variations in the electron identification and add ±0.0007 to 
the combined result. Details of the lifetime fit (binning, fit range, and D° direction choice) 
add ±0.0004. Variations in particle identification selections and details of the refit procedure 
add ±0.0004. Other systematic effects are studied and found to be negligible. 

All the systematic uncertainties are combined in quadrature to yield a total systematic 
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TABLE I: Summary of uncertainties on Rm- 



Source 


Size 


Statistics 


±0.0029 


Simulation Statistics 


±0.0023 


Decay Time Shape 


±0.0014 


Q Parameterization 


±0.0008 


Electron Identification 


±0.0007 


Decay Time Fit Details 


±0.0004 


Particle ID and Refit Details 


±0.0004 


Systematic Total 


±0.0029 



uncertainty of ±0.0029. All of the uncertainties are summarized in Table |T] 

We see no sign of D mixing in either channel and set limits on Rm based on our central 
value of 0.0016 ± 0.0029 ± 0.0029 where the statistical and systematic uncertainties are 
combined quadraticaliy. We assume the uncertainty follows a Gaussian distribution and 
exclude the unphysical region, Rm < 0. This gives Rm < 0.0078 at 90% C.L. and Rm < 
0.0091 at 95% C.L. We agree with previous measurements and set limits as expected given 
the size of our data sample. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation and the U.S. Department of Energy. 
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